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INTRODUCTION 


Based  on  the  investigations  of  Barsoum,^-  Henshell  and  Shaw, 2  quarter- 
point  quadratic  elements  have  been  successfully  used  as  crack  tip  elements  in 
fracture  mechanics.  This  concept  of  singular  element  was  extended  to  cubic 
isoparametric  elements. ^  Recently  it  was  discovered  by  Lynn  and  Ingraffea^ 
that  under  special  configuration,  transitional  elements  improve  the  accuracy 
of  stress  intensity  factor  computations.  These  transitional  elements  are 
located  in  the  immediate  vicinity  of  the  singular  elements  with  the  mid-side 
nodes  adjusted  as  to  reflect  or  extrapolate  the  square  root  singularity  on  the 
stresses  and  strains  at  the  tip  of  the  crack. 

In  this  report,  we  have  obtained  the  locations  of  mid-side  nodes  of  these 
transitional  elements  for  the  quadratic  as  well  as  cubic  elements.  Explicit 
computations  for  a  typical  element  are  symbolically  carried  out  using 
MACSYMA* *.^  These  computations  reveal  that  in  addition  to  the  desired  square 
root  singularities,  the  crack  tip  senses  a  stronger  singularity,  i.e.,  of 
order  one.  Further,  the  strength  of  this  singularity  cannot  be  controlled,  as 
was  possible  for  the  cubic  and  quadratic  collapsed  elements,  where,  by  tying 
the  collapsed  nodes  together,  we  could  easily  abolish  this  strong  singularity. 

These  cubic  elements  also  have  Hibbit-type^  singularities.  The  locations 
of  mid-side  nodes  for  these  singularities  have  also  been  determined. 


References  are  listed  at  the  end  of  this  report. 

*MACSYMA  is  a  large  program  for  symbolic  manipulation  at  HIT. 
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The  cubic  transitional. elements  were  used  for  a  double-edge  crack  prob¬ 
lem,  and  it  was  found  that  there  was  improvement  in  accuracy  for  a  configura- 

*  ■*?'**, 

tion  which  consisted  only  of  singular  and  transitional  elements.  However,  for 
a  well  laid  out  grid,  the  improvement  was  only  marginal.  MACSYMA  has  proved 
to  be  an  indispensable  tool  for  the  present  investigation. 

SECTION  I 

Consider  a  quadratic  quadrilateral  isoparametric  element, 

8  8 

x  -  I  NtXi  ,  y  -  I  NiYi  (1) 

i-1  i-1 

8  8 

u  -  l  NiUi  ,  v  -  l  NiVi  (2) 

i-1  i-1 

where  Nj  are  the  shape  functions  of  'Serendipity'  family, ^  and  are  given  by, 
CORNER  1 

NODES  :  N1  ■  7  (l-5)(l-n)(-S-n-l)  ,  etc.  (3) 

A 

MID-SIDE  1 

NODES  :  N5  “  ”  ( 1— C2) (1— n)  ,  etc.  (4) 

Without  loss  of  generality  consider  the  sectorial  element,  together  with  the 
mapped  unit  element  in  the  transformed  plane,  shown  in  Figure  1.  For 
simplicity,  considering  the  one-dimensional  case  along  line  1-2  in  Figure  1 
(i.e.,  n  -  -1)  we  have  from  (1) 

1  1  . 

x  -  -  U5-1)  +  -  CC1+OL  +  (1-C2)8L  (5) 

2  2 
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TRANSITION  ELEMENT 
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FIGURE  I.  QUADRATIC  QUADRILATERAL  ISOPARAMETRIC 
ELEMENT  AS  TRANSITION  ELEMENT 
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The  condition  for  the  coalescence  of  roots  of  (5)  at  x  ■  0,  together  with  the 
condition  that  PL  >  1  gives 


PL 


L+2/L+1 


(6) 


This  is  the  result,  in  a  slightly  different  form,  obtained  by  Lynn  and 
Ingraffea.^  With  this  location  of  mid-side  nodes,  the  mapping  of  the  general 
element  of  Figure  1  becomes,  from  (1)  and  (2), 

x  -  -  {(n+1)  cos  a  +  (1-n)} {5(/l-1)  +  (/l+1)}2  (7  j 

8 

y  ■  -  (n+l){c(/L-l)  +  (/l+1)>2  sin  a  (8) 

8 

The  Jacobian  of  the  transformation  (1)  and  (2)  is  then  given  by 

3  ( x  y )  1 

J  - - -  —  (/l-1) { £(/L-l)  +  (/L+1)}3  sin  a  (9'i 

3(C,n)  16 

As  can  be  seen  from  (7),  (8),  and  (9),  the  Jacobian  has  a  third  order  zero 
while  x  and  y  have  second  order  zeroes  at 


/L+1 

/L-l 


(io; 


Using  the  inverse  of  the  Jacobian  matrix,  the  strain  component  can  be  written 


as 


3u  1  9u  dy  3u  dy, 

3x  J  3£  dn  3n  d E, 


(ID 
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Substituting  the  various  derivatives  and  collecting  terms  we  get 

9u  A1  A2 

=, - + - +  A3  (12) 

8x  (CC/L-l)  +  /L+1)2  (S(/l-1)  +  /L+l) 

where  A^  ,  A2,  and  A3  are  given  in  Appendix  A. 

Comparing  (12)  with  (7)  and  (8)  it  is  seen  that  the  strain  component  not 
only  has  singularity  of  order  one  half  but  also  of  order  one*  Similarly  we 
have 

3u  1  ^  3u  dx  +  3u  dXj  ^4  ^5  ^  ^ 

3y  J  9 5  dn  3n  d£  (C(/l-1)  +  /l+1)2  (5(/l-1)  +  /l+1)  6 

where  A4,  A5,  Ag  are  given  in  Appendix  A. 


SECTION  II 

Consider  now  the  cubic,  12-node,  quadrilateral  isoparametric  element, 


and  displacements 


12  12 

X  -  1  NiXi  ,  y  -  l  NiYi 

i=*l  1*1 


12  12 

u  “  l  ,  v  -  l  NiVi 

i=*l  i=*l 


(14) 


(15) 


where  the  shape  functions  are  given  by 


CORNER  1  .  .  .  , 

NODES  :  N1  =“  ~~  (l-O(l-n){9(£2+n2)-10}  ,  etc. 


(16) 


MID-SIDE  9 

NODES  :  N2  =  (l-30(l-S2)(l-n)  >  etc.  (17) 
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The  general  transitional  element  together  with  its  map  in  5-n  plane  is  given 
in  Figure  2.  For  simplicity  consider  the  one-dimensional  case  along  line 
1 -2-3-4  (i.e.,  n  -  -1), 


x 


1 

16 


U3 (-9+2701L-27S2L+9L)  +  52(9-901L-9B2L+9L) 


+  ( 1-270 XL+2702L-L)  +  (-1+901L  +  902L-L)) 


(is;. 


The  requirement  that  (18)  be  quadratic  in  K,  together  with  the  condition 
of  coalescence  of  roots  gives  the  following,  physically  possible  solution  for 
locations  of  mid-side  nodes  for  all  L, 


0|L 


L+4/L+4 


4L+4/L+-1 

02L  - - 

9 


(19) 


With  the  above  values  the  general  mapping  of  the  element  shown  in  Figure  2 
then  becomes 

x  -  -  ((n+l)cos  a  -  (n-l))U(/L-l)  +  (/L+l)}2  (20) 

8 

y  -  -  (n+l){^(/L-l)  +  (/u-1)}2  sin  a  (21) 

8 

and  the  Jacobian  of  the  transformation  becomes 

j  -  ti-lll  -  L  (/l-1)(5(/L-1)  +  (/W-l))3  sin  a  (22) 

9(5,n)  16 

These  expressions  are  the  same  as  for  quadratic  elements  (compare  eqs.  (7), 
(8),  and  (9)),  and  hence  the  Jacobian  has  third  order  zeroes  and  x,y  have 
second  order  zeroes,  at 

/iffl 

5 . -  (23) 

/L-l 


6 


SINGULAR  LINE 


(1,-1) 


FIGURE  2.  CUBIC  QUADRILATERAL  ISOPARAMETRIC 
ELEMENT  AS  TRANSITION  ELEMENT 
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Following  the  procedure  outlined  before,  the  strain  components  can  be  obtained 
from  the  following 


3u  B1  ®2 

—  - - + - +  B3  , 

3x  (£(/!-!)  +  /LH)2  (£(/I-l)  +  /U-l) 


3u 


(24) 


(25) 


where  Bj  through  B5  are  given  in  Appendix  A.  Similar  expressions  hold  for 
derivatives  of  v.  Equations  (24)  and  (25)  again  reveal  the  same  kinds  of 
singularities  as  (12)  and  (13). 

SECTION  III 

In  the  cubic  elements  there  is  an  additional  set  of  locations  of  mid- side 
nodes  which  give  Hibbit-type^  singularity.  This  is  obtained  from  the  condi¬ 
tion  that  all  the  three  roots  of  (18)  coalesce.  The  location  of  nodes  is 
given  by 


Ll/3  +  2  3 

piL  -  c~ — > 


(26) 


and  the  transformations  become 


(27) 


8 


and  the  Jacobian  becomes 


(28) 


Following  the  procedure  outlined  before,  it  can  be  shown  t^hat 


3u 


C2 


3x 


(sa1/3-!)  +  L1/  3+l )  3  (e(L1/3-l)  +  L1/^!)2 


C3 


+ 


(CCL1/3-!)  +  lJ3+1  )  4 


(29) 


The  above  equation  indicates  that  in  this  case  the  singularities  are  of  order 
1,  2/3,  and  1/3.  This  combination  is  of  no  immediate  interest  in  linear 
fracture  in  homogeneous  media. 

SECTION  IV 

The  sample  problem  of  a  double-edge  cracked  plate  of  Ref.  4  was  selected 
for  numerical  assessment  of  transition  elements  when  used  with  12-node  col¬ 
lapsed  singular  elements.  Figure  3  is  an  idealization  we  usually  take  for 
such  a  mode  I  crack  problem.  The  distance  p  between  the  crack  tip  and  the 
nearest  node  in  a  collapsed  element  is  often  taken  in  the  range  of  0.5%  to  3% 
of  the  crack  length  a.  The  ratios  a/b  and  b/c  are  usually  in  the  range  of  2 
to  10.  Stress  intensity  factors  for  several  values  of  p,  b/c,  and  a/b  with 
and  without  the  use  of  transition  elements  are  tabulated  in  Table  I.  Compar¬ 
ing  to  the  reference  value,  Kj  »  a/rraF(a/2a)  ,  where  F ( 1  / 2 )  =*  1.184,7  t^e  per¬ 
centage  errors  A%  are  also  shown  in  the  table.  The  result  with  the  use  of 
transition  elements  is  better  only  when  a  very  large  ratio  of  b/c  ( =20)  is 
used. 
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a  =  1.8 


c  =  9p 


FIGURE  3.  AN  IDEALIZATION  FOR  A  QUARTER 

OF  A  DOUBLE -EDGE  CRACKED  PLATE 
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TABLE  I.  STRESS  INTENSITY  FACTOR  AND  PERCENTAGE  ERROR  FOR  A  DOUBLE-EDGE 
CRACKED  PLATE  USING  12-NODE  COLLAPSED  SINGULAR  ELEMENTS  WITH 
AND  WITHOUT  TRANSITION  ELEMENTS.  FINITE  ELEMENT  IDEALIZATION 
OF  FIGURE  3. 


- r 

1  1 

|  Without  Transition 
Elements 

|  With  Transition  [ 
Elements 

1  1 
p 

b/c 

a/b 

1 

|  SIF 

A% 

1  SIF 

A%  | 

i  r 

1  0.005  1 

1  1 

4 

10 

!  2.8808 
i 

2.31 

I  2.8736 

2.06  | 

|  1 

1  1 

10 

4 

1 

|  2.8376 

1 

0.78 

2.7831 

-1.16  1 

1  1 

1  1 

20 

2 

1 

1  2.9863 

1 

6.06 

I  2.7851 

-1.09  | 

1  1 

|  0.01  1 

1  1 

4 

5 

1 

I  2.7986 

-0.61 

|  2.7926 

-0.82  | 

1 

1  1 

1  1 

10 

2 

i 

1  2.8334 

1 

0.63 

1  2.7813 

-1.22  1 

TABLE  II.  STRESS  INTENSITY  FACTOR  AND  PERCENTAGE  ERROR  FOR  A  DOUBLE-EDGE 
CRACKED  PLATE  USING  12-NODE  COLLAPSED  SINGULAR  ELEMENTS  WITH 
AND  WITHOUT  TRANSITION  ELEMENTS.  FINITE  ELEMENT  IDEALIZATION 
OF  FIGURE  4. 


i - r 

1  1 

1  Without  Transition 

Elements 

With  Transition 
Elements 

1  1 

P  ] 

a  /c 

|  SIF 

A% 

|  SIF 

r 

A%  | 

1 

|  0.005  | 

1  | 

40 

|  3.325 

18.09 

I  2.7658 

-1.77  I 

j 

1  1 

1  o.oi  I 

1  1 

20 

2.963 

5.23 

1  2.7654 

-1.79  | 

| 

1 

1  0.02  I 

1  | 

10 

1  2.8115 

-0.15 

1  2.7650 

-1.80  | 

| 

1 

I  0.04  | 

1  1 

5 

|  2.7632 

-1.86 

I  2.655 

-5.71  | 

_ L 

11 


Another  idealization,  Figure  4,  similar  to  the  one  used  by  Lynn  and 
Ingraffea^  was  used  to  recompute  stress  intensity  factors  for  various  values 
of  a/c  to  see  whether  the  transition  elements  in  cubic  isoparametric  elements 
can  give  as  good  improvement  in  accuracy  as  reported  in  Ref.  4  in  the  quadrat¬ 
ic  isoparametric  case.  These  results  are  tabulated  in  Table  II.  It  shows 
again  the  result  obtained  from  the  use  of  transition  elements  is  better  only 
when  a  very  large  ratio  of  a/c  is  used. 

In  this  report  the  stress  intensity  factors  were  calculated  from  the  nor¬ 
mal  component  of  displacement  of  the  node  on  the  crack  surface  and  nearest  to 
the  crack  tip.  It  usually  gives  better  results  than  the  average  value  comput¬ 
ed  from  nodal  displacements  along  the  rays  from  the  crack  tip  at  various 
angles.® 

For  elastic  crack  problems,  the  correct  order  of  singularity  at  the  crack 
tip  is  taken  care  of  by  the  collapsed  singular  elements.  The  use  of  transi¬ 
tion  elements  does  not  practically  improve  the  accuracy. 

CONCLUSIONS 

In  this  report  we  have  been  able  to  obtain  explicit  expressions  for  sin¬ 
gularities  the  crack  tip  senses  from  a  transitional  element.  The  application 
of  these  elements  for  a  few  practical  problems  of  fracture  mechanics  as  well 
as  stress  concentration  factors  has  been  partially  successful.  It  is  believed 
this  is  because  the  crack  tip  senses  not  only  the  square  root  singularity,  but 
also  a  stronger  singularity.  The  strength  of  this  singularity  cannot  be  con¬ 
trolled  as  was  possible  for  collapsed  singular  elements,  where  the  strong  sin¬ 
gularity  was  essentially  eliminated  by  tying  the  nodes  together. 
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FIGURE 


.  A  SIMILAR  IDEALIZATION  USED  IN  [<]  FOR  A 
QUARTER  OF  A  DOUBLE-EDGE  CRACKED  PLATE 
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APPENDIX  A 


In  this  appendix  we  give  the  explicit  expressions  for  the  coefficients  of 
the  various  terms  in  the  strain  components  given  in  the  text. 

R  =  /L 

2(n+l)  , 

Ai  =  ----"2  {4n(R~l) [Rug-ug]  +  4R[U7~U5] 

+  (2nR+R-2n+l) [u2“Ru4]  +  (2nR-R-2n-l) [U3-RU1] } 

1 

A3  ”  -  --”"2  { 2(3n2+4n+l) (R-l) (ug-ug)  +  4(n+l) (R+1)u7 

-  4(n+3)(R+l)u5  +  (R(3n2+7n+4)  -  3n(n+l))u4 
+  (3Rn(n+l)  -  (3n2+7n+4))u3  +  (R(3n2+5n+4)  -  (3n2+n-8))u2 
-  (R(3n2+n-8)  -  (3n2+5n+4))ui } 

2 

A3  "  -  (-”-2  (2u5”u2-u1) 


2((n+l)cos  a  -  (n-l))  , 

Aa  -  - r - { 4Rn (R-l )ug  +  4R(u7~U5)  -  4n(R-l)ug 

(R-l)zsin  a 


-  (R(2n+1)  -  (2n-l))(Ru4-u2)  +  (R(2n-1)  -  (2n+l))(u3-Ru1)} 


Ac  - - — - - { 2(R-1) [cos  a(3n2+4n+l)  -  (3n2-4n+l)] (ug-ug) 

2 (R-l) 2 sin  a 

+  4(R+l)(cos  a(n+l)  -  (n~3))u7  -  4(R+l)(cos  a(n+3)  -  (n-l))u5 
-(cos  a[R(3n2+7n+4)  -  3n(n+l)]  -  [R(3n2-n-8)  -  (3n2-5n+4)] )u4 
+  (cos  a[3Rn(n+l)  -  (3n2+7n+4)]  -  [R(3n2-5n+4)  -  (3n2-n-8)] )u3 
+  (cos  a[R(3n2+5n+4)  -  (3n2+n-8)]  +  [-3Rn(n-l)  +  (3n2-7n+4)] )u3 
-  (cos  a[R(3n2+n-8)  -  (3n2+5n+4)]  +  [-R(3n2-7n+4)  +  3n (n— 1 ) ] )ui} 


A-l 


a6 


(R-l)2sin  a 


{cos  a(-2u5+U2+u^)  +  (2117-114-113)} 


B1  “  ”:”3  {”  [R2 C 2 7ti2— 18ri — 1 )  -  R(54n2-36n-38)  +  27n2-18n-l]  (R14-U4) 

+  9R(2R+l)(u9-u2)  +  9R(R+2) (u3-ug) 

9  9 

+  -  C R—  1 ) 2  (9ri2-2r|-3)  (u5-Ru12)  -  -  (R-l)  2(9n2+2n-3)  (u6-RUl  x) 


1 


+  -  [R2  (27n2+18ri-l)  -  R(54n2+36n-38)  +  27n2+18n-l]  (u7-Ru10) } 


b2  “  7“-73  {-  7-  [R2(45n3+27n2-n+105)  -  R(90n3+54ri2-146ri-222) 

CK— 1)  16 

9 

+  45n3+27ri2-37n-3]u1  +  -  (2R2+6R+1)  [  (n+3)u2  -  (ti+1)u9] 

9  , 

-  -  (R2+6R+2)  [  (ri+3)u3  -  (ri+l)uo]  + 

4 

+  77  [R2(45n3+27n2-37n-3)  -  R(90n3+54ri2-146r|-222)  + 

1 6 

+  45ti3+27ti2-ti+105]u4 

9 

-  --  (R-l)2(n+l)[(15ri2-7)u5  +  ( 15n2+6n-5)Ul x] 


+  --  (R-l)2(n+i)t(i5n2+6n-5)u6  +  (I5n2-7)u12] 


-  —  (n+l)[R2(45n2+36n-l)  -  2R(45n2+36n-37)  +  45n2+36n+35]u7 
16  ' 


n+1 


+  ---  [R2(45n2+36n+35)  -  2R(45n2+36n-37)  +  45n2+36n-l] u10} 


Bo  B  -  — I — -  { (2R+1 )ui  -  (5R+4)uo  +  (4R+5)u3  -  (R+2)u4} 

J  2  (R-l)3 

Ba  -  I------------------  {-  [R2(27n2-18n-l)  -  R(54n2-36n-38)  + 

4  sin  a (R-l) 3  4 

+  27n2-18n-l] (-R14+U4)  +  9R(2R+1)(u2-u9)  + 

9 

+  9R(R+2)  (-u3+ug)  +  -  (R-l)2(9n2-2ri-3)  (-U5+RU12) 

4 

9 

+  -  (R-l)2(9n2+2n-3)(u6-Ru11)  + 

4 

+  -  [R2(27n2+18n-l)  -  R(54n2+36n-38)  +  27n2+18n-l] (-u7+Ru10) } 

4 

1  1 

Be  - - r  { —  [R2(45n3(cos  a-1)  +  27n2(cos  a+3)  -  n(cos  a+71) 

sin  a (R-l)3  16 

+  35(3  cos  a+1))  +  R(90n3(-cos  a+1)  -  54n2(cos  a+3)  + 

+  2n(73  cos  a-1)  +  74(3  cos  a+1))  +  45  n3(cos  a-1)  + 

+  27n2(cos  a+3)  -  n(37  cos  a+35)  -  3  cos  a-l]ui 

9 

-  -  ((n+3)co8  a-n+1 ) [ (2R2+6R+1 )u2  ~  (R2+6R+2)u3l 

4 

-  —  [R2  (45n3  (cos  a-1)  +  27ri2(cos  a+3)  -  n ( 3 7  cos  a-35) 

16 

-  (3  cos  a+1))  +  R(90o3(-cos  a+1)  -  54n2(cos  a+3) 

+  2n(73  cos  a-1)  +  74(3  cos  a+1))  +  45n3(cos  a-1)  + 

+  27n2(cos  a+3)  -  n(cos  a+71)  +  35(3  cos  a+l)]u4 


9 

+  —  (R-l)2 (15n3 (cos  a-1)  +  3n2(5  cos  a+7)  -  n(7  cos  a+1) 
16 


-  7  cos  a-5)ue  -  (15n3(cos  a-1)  +  3n2(7  cos  a+5) 

16 


A- 3 


+  n(cos  a+7)  -  5  cos  a-7)(R-l)2u£  + 

+  --  [R2 (45n3 (cos  a-1)  +  27n2(3  cos  a+1)  +  n(35  cos  a+37) 
1 6 

-  (cos  a+3))  +  R(90n3 (-cos  a+1)  -  54n2(3  cos  a+1) 

+  2n(cos  a-73)  +  74(cos  a+3))  +  45n3(cos  a-1)  + 

+  27f|Z(3  cos  a+1)  +  n(71  cos  a+1)  +  35(cos  a+3)]uy 

9 

+  “  ((n+l)cos  a-n+3) [-(Rz+6R+2)ug  +  (2R2+6R+l)ug]  - 
4 

-  —  [R2  (45n3  (cos  a-1)  +  27t)2(3  cos  a+1)  +  (71  cos  a+1)  + 
16 

+  35(cos  a+3))  +  R(90n3(-cos  a+1)  -  54n2(3  cos  a+1) 

+  2n (cos  a-73)  +  74(cos  a+3))  +  45n3(cos  a-1)  + 

+  27nz(3  cos  a+1)  +  n(35  cos  a+37)  -  (cos  a+3)]u^o  + 


+  —  (R-l)2 [15n3 (cos  a-1)  +. 3n2(7  cos  a+5)  +  n(cos  a+7) 
16 

-  (5  cos  a+7)] uu  - 

9 

-  —  (R-l)2 [15n3 (cos  a-1)  +  3n2(5  cos  a+7)  -  n(7  cos  a+1) 
1 6 

-  (7  cos  a+5) ] m2 


Bt  - - r  {(2R+l)[-cos  aui  +  uiq]  +  (5R+4)[cos  au2  -  ug] 

2  sin  a  (R-l)3 

-  (4R+5)[cos  aug  -  uq]  +  (R+2)[cos  au^  -  U7] } 
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